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Abstract: We study a homogeneous and nearly-isotropic Universe permeated by a homo- 
geneous magnetic field. Together with an isotropic fluid, the homogeneous magnetic field, 
which is the primary source of anisotropy, leads to a plane-symmetric Bianchi I model of 
the Universe. However, when free-streaming relativistic particles are present, they gener- 
ate an anisotropic pressure which counteracts the one from the magnetic field such that 
the Universe becomes isotropized. We show that due to this effect, the CMB temperature 
anisotropy from a homogeneous magnetic field is significantly suppressed if the the neutrino 
masses are smaller than 0.3 eV. 
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1. Introduction 

On very large scales, the observed Universe is well approximated by a homogeneous and 
isotropic Friedmann solution of Einstein's equations. This is best verified by the isotropy of 
the Cosmic Microwave Background (CMB). The small fluctuations observed in the CMB 
temperature are fully accounted for by the standard model of structure formation from 
small initial fluctuations which are generated during an inflationary phase. Nevertheless, 
these small fluctuations are often used to limit other processes or components which may 
be present in the early Universe, like e.g. a primordial magnetic field. 

The generation of the magnetic fields observed in galaxies and clusters [|l| is still 
unclear. It has been shown that phase transitions in the early Universe, even if they do 
generate magnetic fields, have not enough power on large scale to explain the observed 
large scale coherent fields 0. These findings suggest that primordial magnetic fields must 
be correlated over very large scales. 

In this paper, we discuss limits on fields which are coherent over a Hubble scale and 
which we can therefore treat as a homogeneous magnetic field permeating the entire Uni- 
verse. We want to derive limits on a homogeneous field from CMB anisotropies. This 
question has been addressed in the past || and limits on the order of B <, 2 x 10 -9 Gauss 
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have been derived from the CMB anisotropies [||]. A similar limit can also be obtained 
from Faraday rotation ||] . 

We show that the limits from the CMB temperature anisotropy actually are invalid if 
free streaming neutrinos with masses m v < T^ ec are present, where T^ ec denotes the photon 
temperature at decoupling. This is the case if the neutrino masses are not degenerate, i.e. 
the largest measured mass splitting is of the order of the largest mass, hence m u £ 0.04eV. 
The same effect can be obtained from any other massless free streaming particle species, 
like e.g. gravitons, if they contribute sufficiently to the background energy density. This is 
due to the following mechanism which we derive in detal in this paper: In an anisotropic 
Bianchi-I model, free streaming relativistic particles develop an anisotropic stress. If the 
geometric anisotropy is due to a magnetic field, which scales exactly like the anisotropic 
stress of the massless particles, this anisotropic stress cancels the one from the magnetic 
field and the Universe is isotropized. Hence the quadrupole anisotropy of the CMB due to 
the magnetic field is erased. This 'compensation' of the magnetic field anisotropic stress 
by free-streaming neutrinos has also been seen in the study of the effects of stochastic 
magnetic fields on the CMB || [| 10 [ for the large scale modes. In our simple analysis 



the mechanism behind it finally becomes clear. 

The limits from Faraday rotation are not affected by our arguments. 

In the next section we derive the CMB anisotropies in a Bianchi I Universe. In Section || 
we show that relativistic free streaming neutrinos in a Bianchi I model develop anisotropic 
stresses and that these back-react to remove the anisotropy of the Universe if the latter is 
due to a massless mode. In Section || we discuss isotropization due to other massless free 
streaming particles, with special attention to a gravitational wave background. In Section || 
we conclude. 

2. Effects on the CMB from a constant magnetic field in an ideal fluid 
Universe 

We consider a homogeneous magnetic field in z— direction, B = Be z in a Universe filled 
otherwise with an isotropic fluid consisting, e.g. of matter and radiation. The metric of 
such a Universe is of Bianchi type I, 

ds 2 = -dt 2 + a\{t){dx 2 + dy 2 ) + a\(t) dz 2 , (2.1) 

where t is cosmic time. The Einstein equations in cosmic time read 

2 ^l + /4l\ 2 = 87rGp> (2 . 2) 

a\\a± \a±J 
Oil Qj I on d I 

J + — + J— = -8ttGP ± , (2.3) 

a || °_L a || a _L 

2 a± + (a±\ 2 = _ 87rGp {2A) 

The dot denotes the derivative with respect to t. We have introduced the total energy 
density p = PB + Pm + Py+Pv + PA, where ps = B 2 /8tt is the energy density in the magnetic 
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field, and p m , p 7 , p v , p\ are as usual the energy densities of matter (assumed to be baryons 
and cold dark matter), photons, neutrinos, and dark energy (assumed to be a cosmological 
constant), respectively. 

All the above constituents of the Universe, except matter (which is assumed to be 
pressureless) also contribute to the pressure components Pn,P±. The contribution from 
the magnetic field is intrinsically anisotropic and given by 

Pb,± = -Pb,\\ = PB , (2.5) 

as can be read off from the corresponding stress-energy tensor. Note that the magnetic 
field B decays as aj 2 , so that ps scales as aj 4 . 

For later reference we define an 'average' scale factor 

a = a/ 3 a^ 3 , (2.6) 

which is chosen such that it correctly describes the volume expansion. 

Let us also introduce the expansion rates H± = a±/a± and Hn = an /an. The 
anisotropic stress of the homogeneous magnetic field sources anisotropic expansion, which 
can be expressed as the difference of the expansion rates, AH = H± — Hn. We combine 
eqs. (|2.4|) and (|2.3| ) to obtain an evolution equation for AH, 



AH + (2H 1 _ + fljj) AH = 8vrG (P± - P\\) ■ (2.7) 
This pressure difference is actually simply the anisotropic stress. More precisely, 
IV = Tf - P6j , P = T?/3 = {2P± + P|| )/3 , 

IT 1 = U 2 2 = P ± -P = ^(P ± -P ll ) , U 3 3 = P ]l -P = -^(P ± -P ll ) . (2.8) 

At very high temperatures, both photons and neutrinos are tightly coupled to baryons. 
Their pressure is isotropic and thus their contribution to the right-hand-side of (|2.7| ) van- 
ishes. The collision term in Boltzmann's equation tends to isotropize their momentum- 
space distribution. Under these conditions the only source of anisotropic stress is the 
magnetic field. The above equation can then easily be solved to leading order in AH, as 
will be carried out in section ||. 

However, as soon as the neutrinos decouple and start to free-stream, their momentum- 
space distribution will be affected by the anisotropic expansion caused by the magnetic field 
and thus they will develop anisotropic stress. As we will show, the neutrino anisotropic 
stress counteracts the one from the magnetic field. This behavior will be maintained until 
the neutrinos become non-relativistic, then their pressure decays. For the temperature 
anisotropy in the CMB it is relevant whether this happens before or after photon decou- 
pling. This depends, of course, on the neutrino masses. 

We introduce the energy density parameters 

8irGp x (t) _ p^ 
' lx[t) - 3Hi(t) ~ Pc (t) ' 
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corresponding respectively to the magnetic field, matter and radiation etc., such that e.g. 

= B 2 /8irp c , Q m = Pm/pc and f2 7 = p 7 /p c - Here we define the 'average' Hubble 
parameter by 



H l 



a± 
a± 



+ 2 



dj_oi| 
a±a\\ 



With this, eq. ( |2.2j ), implies 



fi 7 + Q v + f2 m + $7a = 1 at all times. 



(2- 



(2.10) 



As an alternative, one could have defined the 'average' Hubble parameter as 

1 

3 



H„ 



a i an 
2— + ^ 
an 



It can easily be verified that the difference between these definitions is of the order of the 
small quantity AH = H± — Hu . More precisely, 



H l 



Hi 



| 2 AH 



1 f AH 

3 \~h7 



(2.11) 



We shall mainly use the definition which yields the constraint ( 2.1 0| ) . 

The scaling of the energy densities corresponding to every species follows from the 
stress energy conservation of every single fluid 

— , Pm = Pm\ — ) , PB = P B 



P-y = Pr 



Pu 



(2.12) 



To obtain the above behavior for radiation, it is important to impose that the fluid is ideal, 
i.e. that pressure is isotropic. This is the case if there are sufficiently many collisions, but 
does not hold for free streaming particles as we shall see in the next section. 

At a fixed initial time one may set a± = an as initial condition. Motivated by obser- 
vations, we assume that the scale factor difference always remains small, 

a± — ay 



5< 1 



(2.13) 



To first order in AH <C H, as long as the magnetic field is the only anisotropic component, 
eq. Q2-7D becomes (see also JTT|| ) 

(2.14) 



AH + 3HAH = 8vrG (P ± - P\\) = 6H 2 Q B . 



In the following we consider both Q B and AH as small quantities and want to calculate 
effects to first order in them. To first order, pb oc a -4 oc p 7 . We can therefore introduce 
the ratio 

r = — = —, (2.15) 

which (to first order) is constant. 

In fig. H we plot the scale factor difference 5q — 5 and AH/H as functions of the 
temperature in a first stage where neutrinos, photons and baryons are all tightly coupled 
and the magnetic field is the only source of anisotropy. 
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Figure 1: Temperature evolution of the scale factor difference So — S and AH/H in units of 
r = fls/Q-y when no free-streaming particle compensates the anisotropy produced by the magnetic 
field anisotropic stress. Here So denotes the scale factor difference 5 today. The evolution of the 
'average' scale factor a is the one of a ACDM Universe. As it is shown in section 3, AH/H is 
constant during the radiation dominated era and S is growing. During the matter dominated era 
AH/H is decaying, AH/H ccl/aocT, and S asymptotes to a constant. 



2.1 Lightlike geodesies in Bianchi I 

Let us now determine the CMB anisotropies in a Bianchi I Universe. We are not interested 
in the usual anisotropies from primordial perturbations, which we disregard in our treat- 
ment, but we concentrate on the effect of the global anisotropy, which to leading order will 
result in a temperature quadrupole. 

We choose the tetrad basis eo = dt, e% = aj 1 ^ for i = 1, 2 and e% = a^ l d^. The dual 
basis of 1-forms is given by 6° = dt, 9 l = a±dx l , for i = 1,2 and 9 3 = a^dx 3 . The first 
structure equation, 

d0 a + u a b A 9 b = , 

yields 

^ = ^°, » = 1,2, and uj 3 o = —0°. (2.16) 
a± ay 

The other non- vanishing connection 1-forms are determined by anti-symmetry, cj a fe = —0Jb a . 
After photon decoupling, the photon 4- momentum p = p a e a satisfies the geodesic equation 

^ + <(e b )pV = 0. (2.17) 

Considering the constraint relation for massless particles p a p a = and setting oTq = p° = 
p = W 2^a=\{p 1 ) 2 i where To is a constant with the dimension of energy (or temperature) 
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that multiplies all the components p a , the above equation is solved by 

a,—,—,— , (2.18) 

a± a± aii / 

where n is a unit vector in the direction of the particle momentum and a is determined by 
the condition p a p a = 0. 

n 1 = sin 9 sin <f> , n 2 = sin 9 cos (j> and n 3 = cos 6 . 

The temperature of photons in such an anisotropic Universe for a comoving observer, 
u = dt, is then given by 



T(t,6) = Vab u a p b = p = T a = T J S -^- + ^-^^-[l + 5co S 2 9 + O(5 2 )} . (2.19) 
We set 



-!- f T(t,9)sm6d9d0= — 



to be the temperature averaged over directions. Note that for 5 = and ao = 1, To is 
simply the CMB temperature at time to. For the temperature fluctuations to first order 
in 5 we obtain 

^ - rM j~ T = ^(3cos 2 - 1) + 0(* 2 ) = y20(n) + ° (52) • (2 - 20) 

Hence, to lowest order in 5 a homogeneous magnetic field generates a quadrupole which is 
given by 

-i 2 X67T 

C2 = 5 l° 2 -| 2 = 5 l«2o| 2 = ~ 0.22 x 5 2 . (2.21) 

m=-2 ' ' 

Of course, in principle one can set S(t\) = at any given moment t\ which then leads 
to 4?r-(ti) = 0. However, for the CMB we know that photons start free-streaming only 
at tdec when they decouple from electrons. Before that, scattering isotropizes the photon 
distribution and no quadrupole can develop 1 . In other words, we have to make sure that 
the anisotropy-induced quadrupole is fixed to zero at decoupling and only appears as a 
result of differential expansion between last scattering and today. This can be taken into 
account by simply choosing the initial condition <5(idec) = 0. Without this initial condition 
we have to replace S(t) by 5(t) — 5(td ec ) in eq. ( ^21| ) 2 . The general result for the CMB 
quadrupole today is therefore 

C2 = ^[6(t )-5(t dcc )} 2 . (2.22) 



1 This is not strictly true and neglects the slight anisotropy of non-relativistic Thomson scattering. 

2 More generally, one can say that 5 itself is not a quantity with a physical meaning as long as no reference 
value is specified. In physical terms, only the difference of <5 between two instants of time can be a relevant 
quantity. 
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2.2 The Liouville equation 

At this stage it is straightforward to check that the exact expression found above for the 
temperature, eq. ( 2.19| ), satisfies the Liouville equation for photons (see, e.g. |l2|]) 

p*e a (f 7 )-u i b (p)p b ^L = Q, (2.23) 

when we make the following Ansatz for the distribution function of massless bosonic par- 
ticles in our Bianchi I Universe 



Pi 



= V(P r ) 2 + (P 2 ) 2 , P\\=P 3 , P=jpl+Pj=p\ (2-24) 



iV 7 I 



/ ^ T)= (2vr)3 e ^-l ' T = T ^ 9 )- ( 2 - 25 ) 
Indeed, using eqs. ( [2.16D , we find the following differential equation for the temperature T 

~Jyr~ET l P± ~ — 7^P\\ = • (2-26) 

ol at a± p± an opu " 

With (HH) this can be written as 

r + ±(E±)\?l(n) 2 = . (2 . 27) 

T a±\p ) ay \ p J 

The time behavior of the different components of the photon momentum are given by 
eq. ( 2.1 8|) and one immediately sees that expression ( [2.19 ) for the temperature solves the 
above differential equation. 

Moreover, defining the time dependent unit vectors p l = p % jp and the shear tensor 

°ab = $ab ~ \$ c c Kb , where -& ab = ^ (V a u b + V h u a ) and h ab = rj ab + u a u b , 
one can rewrite the above Liouville equation as 

(p) = -pcTijff , (2.28) 

where p denotes the redshift-corrected photon energy defined as p = ap. This last expres- 
sion agrees with the corresponding equation given in [13]. 

Using the expression for the distribution function of massless fermions, we can also com- 
pute the pressure of neutrinos once they start free-streaming. Indeed, given the fact that 
neutrinos can be considered massless before they become non-relativistic, their geodesic 
equation has the same solution as the one for photons found above, therefore we immedi- 
ately obtain the time behavior of their temperature in an anisotropic Bianchi I background. 
Taking also into account the fact that neutrinos are fermions, their distribution function 
reads 

/VI T 
Ut,T)= » j 3ep/r + 1 , with T (t,0) = ^-[l + 5cos 2 + O(5 2 )] . (2.29) 



-7- 



Note that the parameter T appearing in the neutrino distribution function in not a temper- 
ature in the thermodynamical sense as the neutrinos are not in thermal equilibrium. It is 
simply a parameter in the distribution function and its time evolution has been determined 
by requiring the neutrinos to move along geodesies i.e. to free-stream. 

This distribution function remains valid also in the case where neutrinos are massive, 
i.e. T v < m v . The only difference is that the relation p° = p changes to p° = \fp 2 + (m v a) 2 
which of course affects the momentum integrals for the neutrino energy density and pres- 
sure. 

The energy T u /ciq is the present neutrino 'temperature' in the absence of a homoge- 
neous magnetic field (5 = 0). The energy density p u and the pressure P Uj i in direction i 
with respect to our orthonormal basis are 

Pu = N u J d 3 p f„(t,T)p° (2.30) 
Pu,i = N u j ' d z pf u (t,T)^. (2.31) 



Calculating the integral ( 2.31| ) for relativistic neutrinos to first order in 5 in the direc- 



tions perpendicular and parallel to the magnetic field direction, one finds for the neutrino 
anisotropic stress in the ultra-relativistic limit 

Pu,±-P„,\\ *~Pu(8-8J , (2.32) 
" 15 

where <5* is the value of 5 at neutrino decoupling and can be fixed to zero for convenience. 

The temperature dependence of the neutrino pressure is shown in fig. To leading 
order, this also gives the temperature dependence of the neutrino anisotropic stress. From 
the plot it is clear how the pressure scales as a -4 as long as the neutrinos are ultra- 
relativistic. Once they have become effectively non-relativistic, their pressure decays more 
rapidly, as a -5 . The break in the power law is not precisely at T = m u , but at a somewhat 
lower temperature. Because the neutrinos still have the highly relativistic Fermi-Dirac 
distribution from the time of their thermal freeze-out, it takes some additional redshift 
until they behave effectively non-relativisic. This will have some effect on the estimates 
for the residual CMB quadrupole, as we shall see in sec. ^, in particular the discussion of 
fig- I 

3. Neutrino free-streaming and isotropization 
3.1 Massless free- streaming neutrinos 

We now calculate the effect of free-streaming neutrinos perturbatively, i.e. to first order 



in 5, AH/H and CIb- We linearize eq. (2.7), taking into account the contribution of a 



free-streaming relativistic component to the right-hand side. We have shown that this 
contribution, to leading order in 5, is given by eq. ( [2.32|) . Furthermore, up to 0(5 2 ) 
corrections, 5 is just the integral of AH, 

AH(i0dt' = m^$-ln^4-S-<5* , (3.1) 

t. a l| W a l| (**) 
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Figure 2: Temperature evolution of the neutrino pressure P u .±_ normalized to the neutrino energy 
density p v . The temperature is given in units of the neutrino mass. Note that the break in the power 
law is not at T — m, but at somewhat lower temperature. This is due to the highly relativistic 
Fermi-Dirac distribution of the neutrinos, see also the discussion of fig. 5 in sec. 3.3. 



so that to first order we can identify AH ~ 5. 

Inserting this back into eq. fl2.7|) we find, to linear order in 5, 

) = 6H 2 n B . 



5 + 3H5 + -H 2 n u (8 
5 



(3.2) 



Note that, because we are working at linear order, it is not important with respect to which 
scale factor H, £l u and O b are defined in ( |3.2[) . We will now give analytic solutions to this 
equation for different regimes in the evolution of the Universe. 

Let us begin at very high temperature where the neutrinos are still strongly coupled 



to baryons. In this case they do not contribute to eq. (3.2) since their pressure is isotropic 
(P v> ± — Pj/j ~ 0) given the high rate of collisions. Furthermore, since we are in the 
radiation dominated era (a oc t 1 ^ 2 ), we have H = l/2t, and Qb is constant. The 
solution to eq. (|3.2|) in this case is 



6(t) = AH(t) = *k + _g_ 



(3.3) 



t ' i 3 / 2 ' 

The dimensionless quantity AH/H hence asymptotes to a constant, since the homogeneous 
piece decays like a -1 : 

AH 

m B ■ (3.4) 



H 

AH soon becomes insensitive to the initial conditions and only depends on Qb- This 
also shows that in the absence of an anisotropic source (Ob = 0), the expanding Universe 
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isotropizes. Integrating this equation and remembering that £Ib = constant to first order 
in a radiation dominated Universe, we obtain 



S(t) -5(t') = 3n B Ht/t')- (3.5) 

As the Universe reaches a temperature of roughly 1.4 MeV, the neutrinos decouple and 
begin to free-stream, giving rise to the corresponding term in eq. (|3.2| ). In the radiation 
dominated era, Sl u remains constant as long as neutrinos are ultra-relativistic 3 . This is 
certainly true for temperatures well above a few eV. In this regime, the general solution of 
eq. ( |3,2D is given by 

S( t ) ~$* = ^jf- + t~ 1/4 (c+tV2*V5-l/16 + c _ t -iV*V5=l7l6) (3 6) 

For Q v > 5/32, the homogeneous part is oscillating with a damping envelope oc i -1 / 4 oc 
a -i/ 2 . This means that AH = 5 will decay within a few Hubble times, which is a mat- 
ter of seconds at the temperatures we are talking about. After that, 5 — 5* will remain 
constant at the value of (15/4) Qb/^u until the neutrinos become non-relativistic. Then 
their pressure drops dramatically and so does their anisotropic stress. Until this time, 
the Universe expands isotropically, because the anisotropic stress of the magnetic field is 
precisely cancelled by the one of the neutrinos. Remember that a constant 5 can always 
be absorbed in a re-scaling of the coordinates and has no physical effect. Fig. [3] shows the 
temperature evolution of 5 — 5* in the radiation dominated era from neutrino decoupling 
until T = lOOeV. 

This mechanism rests on two important facts. Firstly, as long as neutrinos are ultra- 
relativistic, they redshift in the same way as the magnetic field, meaning that VLb/£I v is 
constant. Once the anisotropic stress of the neutrinos has adjusted to the magnetic field, 
their sum remains zero independent of the expansion of the Universe which is now in a 
Friedmann phase. Secondly, the efficiency of the effect hinges on the absolute value of £l u . 
In the radiation dominated era (after positron annihilation), we have £l u ~ 0.4 so that 
Q, v > 5/32, and hence the system behaves as an underdamped oscillator with a damping 
envelope oc t -1 / 4 . Had the density parameter of the free-streaming particles been less 
than 5/32, the behavior would be that of an overdamped oscillator. As it is evident from 
eq. Q3.6| ), for Q, v <C 5/32 there would be a mode which decays extremely slowly, roughly 
as £ -4f W 5 . This is why a strongly subdominant free-streaming component cannot damp 
the anisotropy efficiently. As we shall discuss in section |4|, a primordial gravitational wave 
background could play the role of such a free-streaming component if SIgw ^ 5/32. 

3.2 Massive neutrinos 

The neutrinos become non-relativistic roughly at the time when their temperature drops 



below their mass scale. Current bounds on the neutrino mass 14] are such that the 



3 Actually, f2„ changes slightly when electron-positron annihilation takes place, a process which heats up 
the photons but not the neutrinos. This happens at a temperature close to the electron mass. After that, 
~f remains constant until the neutrinos become non-relativistic. 
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Figure 3: Temperature evolution of S — <5» from neutrino decoupling to T — lOOcV. After decou- 
pling, (5 — 5* begins to oscillate with a decreasing amplitude around the constant jf~ > as predicted 



by the analytic solution (3J3). This qualitative behavior is independent of the initial conditions. 



highest-mass eigenstate is somewhere between ~ 1 eV and ~ 0.04 eV. Since the neutrino 
mass splitting is much below 1 eV, an eigenstate close to the upper bound would mean that 
the neutrinos are almost degenerate and hence become non-relativistic all at the same time. 
If this happens before photon decoupling, i.e., if m u > 0.3 eV, the isotropization effect will 
not be present and the CMB will be affected by the anisotropic expansion sourced by the 
magnetic field. However, if the neutrinos remain ultra-relativistic until long after photon 
decoupling, the CMB quadrupole due to anisotropic expansion will be reduced because the 
neutrinos maintain expansion isotropic until they become non-relativistic. 

In order to quantify this statement, we repeat the above calculations for the matter 
dominated era. For our purposes, this is a reasonable approximation for the time between 
photon decoupling and today. At decoupling, radiation is already subdominant, and on 
the other hand vacuum energy only begins to dominate at redshift z ~ 0.5. We therefore 
expect that both give small corrections only. 

For completeness, we also give the solution of eq. ( [3.2| ) in a matter dominated Universe 
for the case where we ignore any contributions from free-streaming particles (neutrinos 
and, after decoupling, also photons). During matter domination we have H = 2/3t and 
£Ib oc a -1 oc i~ 2 / 3 . The solution to ( p. 14 ) hence reads 

8Q B (t) C 
t + t 2 



5(t) = AH(t) = ^ + J . (3.7) 



The homogeneous mode again decays more rapidly than the particular solution, so that 
the dimensionless quantity AH/H is again asymptotically proportional to £Ib- Instead of 



eqs. fl3^4|) , (p\q) , we have 

A# '* 



if 



12ft B , 5(t) - S(t cq ) = j AHdt ~ 12 [0 B (teq) - fi B (*)] • ( 3 - 8 ) 



Let us now take into account a free-streaming component. We want to estimate the 
effect on the photon distribution function caused by anisotropic expansion in two cases. 
Case A: the neutrinos become non-relativistic before photon decoupling. Case B: the neu- 
trinos become non-relativistic after photon decoupling. As an approximation, we assume 
that this happens instantaneously to all neutrino species, such that the contribution of 
neutrinos to eq. ( |3.2| ) disappears abruptly. We know that the neutrinos are in fact spread 
out in momentum space and also have a certain spread in the mass spectrum, so in reality 
this will be a gentle transition. However, we only want to estimate the order of magnitude 
of the effect and are not interested in these details at this point. More precice numerical 



results will be presented in sec. 3.3. Let us consider case A first. 



3.2.1 Case A: neutrinos become non-relativistic before photon decoupling 

We know that AH is very nearly zero when the neutrinos become non-relativistic. After 
that, AH/H will start to grow again to approach the value 6£Ib during radiation domina- 
tion and 12Qb during matter domination. As boundary condition at photon decoupling, 
we will hence assume AH/H = xQb with x < 12. This number can in principle be 
computed given the neutrino masses and the evolution of the scale factor across matter- 



radiation equality. We shall solve the full equations in subsection 3.3; here we just want to 
understand the results which we obtain there by numerical integration. The free-streaming 
component we are interested in now are the photons after decoupling. We therefore iden- 
tify <5* = 5(£dec)j where tdec denotes the instant of photon decoupling. Furthermore, in 
eq. ( |3.2j ) we replace Q u by f2 7 , our new free-streaming species. With f2 7 oc t~ 2 / 3 in the 



matter dominated era, the (not so obvious) analytic solution to eq. (|3.2| ) is 

S(t) - 5(t dcc ) = + C [f(t) cos f(t) - sin f(t)] + D [f(t) sin f(t) + cos f(t)] , (3.9) 

where we have introduced f(t) = 4y / 2r2 7 (i)/5. The time derivative of eq. fislj ) yields 

^ = ^ 7 [C sin f(t) - D cos f(t)] . (3.10) 

Note that the slowly decaying mode has the same asymptotic behavior as (3^8) - in the 
matter dominated era, the free-streaming radiation can never catch up to the magnetic field, 
since both fade away too quickly. In other words, this means that free-streaming photons 
are never able to counteract the magnetic field anisotropy in order to isotropize again 
the Universe, even if they represent the main contribution to the background radiation 
energy density, and the reason for this is that they decouple only after the end of radiation 
dominantion. 

In order to estimate the value of 5 today (to), we can simply take the limit of small 
7 (io) <C 1 of (jO)|), Correction terms are suppressed at least by yO^o) ~ 10~ 2 . We find 

15 $7 r 

S(t ) - 5(t dec ) ~ -j-^f + D . (3.11) 
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The constant D is fixed by the boundary conditions at decoupling, given by AH/H = x£Ib 
and 5 = <5(idec)- These boundary conditions translate to 



D 



n B (t, 



dec) 



^7 (^dec) 
^B^dec) 



sin f(t, 



dec) 



/(*dec) 



16' 



15 
T 



16 



X cos f(t, 



dec J 



^7 (^dec I 



15 / 2x\ 
"T + ( 4 + T j ^ dec) + (^7^ec)) 



(3.12) 



In order to obtain the essential behavior we have expanded the boundary term as a Taylor 
series in 7 (tdec) *C 1> Our final result is 



<5(*o) - <5(t 



dec J 



(3.13) 



up to corrections suppressed by powers of fi 7 (idec)- 

In this case, the CMB quadrupole is not affected by the presence of free-streaming 
neutrinos and we obtain the same result as when neglecting their presence, 



C, 



167T ... . 

225 [<5(Co) 



<K*dec)]' 



768vr 2 
^ ''B^decJ 



O.lV 



(3.14) 



3.2.2 Case B: neutrinos become non-relativistic after photon decoupling 

In this case, the presence of the neutrino anisotropic stress will delay the onset of anisotropic 
expansion until a time t m when the neutrinos become effectively non-relativistic. As before, 
we will ignore that this is a gradual process and simply assume that one can define some kind 
of "effective" t m at which the neutrino anisotropic stress drops to zero. The full numerical 



result is given in section |3.3| . The effect of anisotropic expansion on the photon distribution 
function is estimated as follows. We assume there is no anisotropic expansion between 
photon decoupling and t m . At later times, neutrino anisotropic stress can be ignored. The 
relevant solution ( |3.9| ) is hence obtained with boundary condition 5(t m ) = 0. Working 
through the steps above once again or simply taking the result ( 3.13j ) with idee - > t m and 
x — > 0, one finds 

5(t ) - 5(t dcc ) = 5{t ) - 5{t m ) ~ An B (t m ) . (3.15) 

Since Q,b decays as a -1 , the effect of anisotropic expansion in case B is suppressed by 
roughly a factor of a(tdcc)/(3a(t m )) with respect to case A. For light neutrinos with a 
highest-mass eigenstate close to the current lower bound, this factor can be as small as 
~ 0.03, loosening the constraint on a constant magnetic field from the CMB temperature 
anisotropy correspondingly. Constraints coming from Faraday rotation are not affected. 

Clearly, the heaviest neutrino becomes massive at redshift z m = m u /T u > 0.04eV/T Jy ~ 
200. One might wonder whether isotropization can be supported even if only one neutrino 
remains massless, since its contribution to the energy density is Vt u \ ~ 0.23S1 7 . The problem 
is however that, as soon as one neutrino species becomes massive, the equilibrium between 
the magnetic field and the neutrino anisotropic stresses is destroyed and, as we have seen 
under case A, where one still has free streaming photons, it cannot be fully re-established 
in a matter dominated Universe. 
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3.3 Numerical solutions 

In order to go beyond the estimates derived so far, we have solved eqs. (pT2 -2.4) numerically 
with cosmological parameters corresponding to the current best-fit ACDM model [15]. We 
use cosmological parameters = 0.73, Q m = 0.27 today, where f2 m includes a contribution 
of massive neutrinos 4 which we approximate by £l u h 2 = N„m u /94:eV with N v ~ 3. The 
contribution to the right-hand side of eq. (2.7) from free-streaming neutrinos is obtained 
by integrating eq. ( [2.31 ) with the full distribution function for massive fermions. More 
precisely, we compute the full distribution function to first order in 5 and perform the 
integration numerically, including the neutrino mass as a parameter. We begin to integrate 
deep inside the radiation dominated era, when the neutrinos are still relativistic but already 
free-streaming. The asymptotic behavior of solution can be used as initial condition at 
neutrino decoupling. The constraint equation (2.2) provides the remaining initial condition. 
We then integrate until the desired time. We define today to by a(to) = 1. 

In fig. |3|, we present the results of the numerical integration from neutrino decoupling 
until today. We plot both 5 — 5* and AH/H in units of the parameter r = fi#/J) 7 so that 
the plots are valid for arbitrary magnetic field strengths, as long as r <C 1. After neutrino 
decoupling, 5 oscillates and reaches its constant value as in eqs. ( |3.6| ), (|3.9|), while AH = 5 
oscillates and decays. We choose as initial condition 5 = 5* = at neutrino decoupling. 
Once the temperature of the Universe reaches the neutrinos mass scale, neutrino pressure 
decreases and they become non-relativistic. At this point, they can no longer compensate 
the anisotropic pressure of the magnetic field, and both 5 and AH begin to grow. How- 
ever, it is clear from fig. Q (upper plot) how, once neutrinos become non relativistic after 
photon decoupling (case B), the growth of 5 is suppressed with respect to case A, where 
this happens before photon decoupling. Moreover, the solid black line in the lower plot 
represents the temperature evolution of AH/H in the case where only the magnetic field 
sources the anisotropy: this makes clear how the absence of any free-streaming particle 
able to counteract the magnetic anisotropic stress leaves the anisotropy of the Universe 
free to grow with respect to its value today. 

Our quantitative final result is shown in fig. p|, where we plot the value of the quadrupole 
generated by a constant magnetic field, rescaled by r 2 , as function of the neutrino mass. 
We weight the final C2 with respect to the quadrupole obtained without considering the 
isotropization induced by free-streaming particles, in order to underline the relative impor- 
tance of this effect. These results clearly show that the CMB quadrupole is significantly 
reduced by neutrino free-streaming only if their mass is smaller than the temperature 
at photon decoupling, m u < T^ cc ~ 0.26 eV. In fact, for neutrino masses in the range 
0.3 eV < m u < 3eV, the quadrupole C2 is reduced by less than a factor 100 from the result 
without a free-streaming component, whereas for < m„ < 0.3 eV, it decreases by several 
orders of magnitude. Note, however, that the effect is not negligible even in the former 
case with relatively large neutrino masses. Fig. || also shows our analytical estimation for 

4 CMB observations actually constrain the matter density at decoupling, such that neutrinos with m„ < 
0.3eV, which are still relativistic at that time, do not contribute to the measurement of f2 m . However, since 
their density parameter today is then also very small, their contribution to the matter density remains 
practically irrelevant. 
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Figure 4: Temperature evolution of AH/H and 6 — <5* for different neutrino masses. We chose 
the initial conditions to be given by <5* = at neutrino decoupling. The black solid line in the 
lower plot represents the temperature evolution of AH/ H in the case where only the magnetic field 
sources the anisotropy and no free-streaming particle is present to compensate this effect. The 
dotted vertical line indicates the instant of photon decoupling. 



the final amplitude of the CMB quadrupole produced by this effect as given by eq. ( 3. 15| ) . 
Of course the value of eq. ( |3.15 ) depends on the time at which neutrinos become effectively 




0.01 eV 



1.0 eV 



Figure 5: Effect of free-streaming neutrinos with different masses on the quadrupole generated 
by a homogeneous magnetic field, weighted on the quadrupole obtained without considering the 
effect of any free-streaming particles. The solid black line represents the result of the numerical 
integration, the dashed blue and red lines correspond to our analytical prediction given by eq. (3.15) 
for two different choices of t m , the time at which neutrinos are effectively non-relativistic (see the 
text for clarification). 



non-relativistic, t m . Once we choose t m to be given by the time at which T = m u , we 
overestimate the final quadrupole amplitude still by one order of magnitude (dashed blue 
line). This is a consequence of the fact that the neutrino distribution function is highly 
relativistic and therefore it takes a further redshift for them to start behaving effectively 
as massive pressureless particles. This has been considered in the more elaborate estimate 
given by the dashed red line where we fix the time t m to be given by the time at which 
d 3 P v I 'fi(lnT) 3 = 0, i.e. the time at which the pressure reaches the break in the power law. 
This is in excellent agreement with the numerical results. 



4. A gravitational wave background and other massless free-streaming 
components in an anisotropic Universe 

From our previous discussion it is evident that any massless free-streaming particle species 
X can isotropize the Bianchi I model with a constant magnetic field, if present with suf- 
ficient contribution 0.x already in the radiation dominated era. This has to be accounted 
for if we want to estimate the CMB quadrupole induced by a homogeneous magnetic field. 

So far we have discussed the standard model neutrinos as an example of such a par- 
ticle. However, also other massless particles can play this role, for instance gravitons, but 
also particle species outside of the spectrum of the standard model. Interestingly, the cur- 



rent bounds on the number of relativistic degrees of freedom during nucleosynthesis, often 
parameterized by the effective number of additional neutrino species AN„, allow for the 
possibility that such a species could be sufficiently abundant. The present bound on N v 
from nucleosynthesis is [14] 



N v = 3.2 ± 1.2 , 

7N ( 4 \ 4//3 
g. = 2+-fl-\ = 3.36 + (N u - 3) x 0.454 

= 3.36 + (0.2 ± 1.2) x 0.454 at 95% confidence. (4.1) 

Here we have taken into account that the photon and neutrino temperatures are related by 
T v = (4/ll) 1//3 r 7 [16]. The effective g* from 7 and three species of neutrino corresponds 



to <7*(7, 3f) = 3.36. This is equivalent to a limit on an additional relativistic contribution 
at nucleosynthesis of tlx ^ 0.2. Prom the solution ( |3 . 6[ ) we know that a free-streaming 
relativistic species with a density parameter tlx ^ 5/32 ~ 0.156 during the radiation 
dominated era will isotropize expansion within a few Hubble times. Since this species will 
presumably decouple before the neutrinos (otherwise it should have been discovered in 
laboratory experiments), expansion can be isotropic already at neutrino decoupling, and 
thus neither the cosmic neutrino background nor the CMB will be affected by anisotropic 
expansion. In this case therefore, unless we are able to detect the background of the species 
X, we will never find a trace of the anisotropic stress produced by a homogeneous magnetic 
field. An interesting example are gravitons, which we now want to discuss. 

Inflationary models generically predict a background of cosmological gravitational 
waves which are produced from quantum fluctuations during the inflationary phase. The 
amplitude of this background, usually expressed by the so-called tensor-to-scalar ratio, tt, 
has not yet been measured, but for a certain class of inflationary models, forthcoming 
experiments such as Planck might be able to detect these gravitational waves. This is in 
contrast to the cosmic neutrino background, for which there is no hope of direct detection 
with current or foreseeable technology. However, this background typically contributes 
only a very small energy density, 

^GW,inf/^ 7 - 10- 10 r T , n T <0. 

Only non-standard inflationary models which allow for > can contribute a significant 
background, see []i~8f , 

Gravitational waves can also be produced during phase transitions in the early Uni- 
verse ]jl9f| , after the end of inflation. Such gravitational wave backgrounds can easily 
contribute the required energy density. Let us therefore concentrate on this possibility. 

If the highest energy scales of our Universe remain some orders of magnitude below the 
Planck scale, gravitational waves are never in thermal equilibrium and can be considered 
as free-streaming radiation throughout the entire history. Therefore, if the gravitational 
wave background was statistically isotropic at some very early time, then any amount 
of anisotropic expansion taking place between this initial time and today will affect the 
gravitons in a similar fashion as any other free-streaming component, and therefore our 
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present gravitational wave background would be anisotropic. Loosely speaking, the in- 
tensity of gravitational waves would be larger in those directions which have experienced 
less expansion in total since the initial time when the gravitational wave background was 
isotropic. 

As we have specified above, with the current limits on A.N V , the density parameter of 
gravitons Oqw during nucleosynthesis can be as large as ~ 0.2. At higher temperatures 
(that is, at earlier times), the number of relativistic degrees of freedom increases (more 
particle species are effectively massless), such that Oqw at earlier time can even be larger 5 . 
It is therefore conceivable that gravitons acquire sufficient anisotropic stress to compensate 
the magnetic field and hence take over the role which neutrinos have played in section |3[ 
As already pointed out, in this case, neither neutrinos nor photons will ever experience any 
significant anisotropic expansion, since the Universe remains in a Friedmann phase after 
the gravitons have adjusted to the magnetic field. Of course, gravitons remain relativistic 
for all times and the mass effect which we discussed for the neutrinos does not occur. 

In order to rule out this scenario, it would be very interesting not only to measure 
the background of cosmological gravitational waves but also to determine whether or not 
it shows a quadrupole anisotropy compatible with such a compensating anisotropic stress. 
Or in other words: just as the smallness of the CMB quadrupole is a direct indication 
for isotropic expansion between decoupling of photons and today, the smallness of the 
quadrupole of a gravitational wave background would inform us about the isotropy of 
expansion between today and a much earlier epoch where this background was generated. 



5. Conclusions 



In this paper we have studied a magnetic field coherent over very large scales so that it 
can be considered homogeneous. We have shown that in the radiation dominated era the 
well known Bianchi I solution for this geometry is isotropized if a free streaming relativis- 
tic component is present and contributes sufficiently to the energy density, Qx ^ 5/32. 
This is in tune with the numerical finding |?], ^, that the neutrino anisotropic stresses 
'compensate' large scale magnetic field stresses. A perturbative explanation of this effect is 
attempted in || . Here we explain the effect for the simple case of a homogeneous magnetic 
field: free streaming of relativistic particles leads to larger redshift, hence smaller pressure 
in the directions orthogonal to the field lines where the magnetic field pressure is positive 
and to smaller redshift, hence larger pressure in the direction parallel to the magnetic field, 
where the magnetic field pressure is negative. To first order in the difference of the scale 
factors this effect leads to a build up of anisotropic stress in the free streaming component 
until it exactly cancels the magnetic field anisotropic stress. This is possible since both 
these anisotropic stresses scale like a -4 . 



5 During a transition from gi relativistic degrees of freedom to (72 < fli, the temperature changes from 
T\ to T2. Since entropy is conserved during the transition we have giTf = (72T2. Hence p2 = giT% = 



a-2 



1/3 
ax) T] 

92 



1/3 



/ \ 1/3 

(§2) Pi > Pi- In other words, the energy density of all species which are still in 



thermal equilibrium increases if one reduces the number of degrees of freedom at constant entropy. 
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In standard cosmology this free-streaming component is given by neutrinos. However, 
as soon as neutrinos become massive, their pressure, P v oc a~ 5 , decays much faster than 
their energy density, p u oc a -3 , and the effect of compensation is lost. If this happens 
significantly after decoupling, there is still a partial cancellation, but if it happens be- 
fore decoupling, the neutrinos no longer compensate the magnetic field anisotropic stress. 
Furthermore, a component which starts to free-stream only in the matter era (like e.g. 
the photons) does not significantly reduce the anisotropic stress. Actually, inserting the 
dominant part of the constant D from eq. ( 3.12 ) in ( |3,10 ) one finds 



— = 12fi B , (5-1) 

like without a free-streaming component. 

This cancellation of anisotropic stresses does not affect Faraday rotation. A constant 
magnetic field with amplitude Bq> 10 -9 Gauss can therefore be discovered either by the 
Faraday rotation it induces in the CMB ||, or, if a sufficiently intense gravitational wave 
background exists, by the quadrupole (anisotropic stress) it generates in it. 

Finally, Planck and certainly future large scale structure surveys like Euclid will most 
probably determine the absolute neutrino mass scale. Once this is known, we can infer 
exactly by how much the CMB quadrupole from a constant magnetic field is reduced by 
their presence. 
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